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摘 要 BEARES RUE T Aa fe E RR 控制 十 分 重要 。 在 对 直径 60 mm 的 40CrNiMo KAA 
棒 的 研究 中 ， 有 限 元 模拟 的 淳 火 应 力 及 其 分 布 与 蕊 射线 衍射 测量 的 结果 很 好 的 符合 ， 由 此 论证 了 所 运用 的 慢 度 场 、 组 织 场 和 应 力 场 相 互 耦合 的 数学 
模型 [包括 我 们 建立 的 相 变 塑 ) 的 正确 性 [OM HEMIA TAKI EH pd ne AMARA, dos Y A EL EUIS TER IRL AERE IL 734 
布 规律 及 其 起 因 以 及 不 同 淳 火 介 质 对 淳 火 应 力 的 影响 规律 ， 为 工程 中 工件 淳 火 开 裂 的 分 析 提 供 了 理论 指导 。 
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ABSTRACT Quenching is one of the most important SM processes for improving the mechanical 
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properties of steel components in manufac industry. Whe quenching stress is a source of cracking, which is 
frequently detrimental to steel properties. Therefore, the investigation of quenching stress is very import. r the 
control of distortion, cracking and residual stress distributions of components. In the study of quenching stfess, the 
measurement of stress distribution is necessary to the stress analysis and design of NV process ever, 
in most cases, the cracking of a quenched component is caused by transient stress diirihg quenching, while 
experiment can only measures the final internal stress (residual stress), ratherthan transi tress. As a result, the 
measurement of residual stress associated with finite element simulation (F as been ea direction in 
investigation of quenching stress. In this work, a full through-hardened 40CzNiMo cylinder with 60 mm diameter 
was water-quenched, and cooling curves at three positions a f the radiu cylinder were first measured. Then, 
an optimized heat transfer coefficient as a function of surface temperature was obtained by fitting with the 


measured cooling curves using the trial and error method. Based on an exponent-modified (Ex-Modified) 


资助 项 目 国家 自然 科学 基金 项 目 No. 51371117 
收 稿 日 期 2016-11-14， 定 稿 日 期 2017-01-18 
作者 简介 刘 玉 ， 男 ，1986 年 生 ， 博 士 生 
通讯 作者 陈 乃 录 ，nlchen@sjtu.edu.cn， 主 要 从 事 钢 铁 材料 淳 火 冷 却 技术 与 先进 高 强度 钢 的 研究 
DOI 10.11900/0412.1961.2016.00509 


normalized function describing transformation plasticity kinetics proposed, the thermo-elasto-plastic constitutive 
equations were deduced. The commercial finite element software, Abaqus/Standard, was used to solve the coupled 
temperature field, microstructure field and stress (strain) field. The results indicate that the quenching stress and 
its distribution predicted by FES is well consistent with those measured by X-ray diffraction, which verified that 
models employed in coupling of thermal field, phase transformation field and stress field including transformation 


compressive stress exists in the core and surface of der and the maximum tensile stress exists at subsurface. 


plasticity function proposed are correct. prie wer of residual stress distribution were revealed that 
The separated calculation of thermal stress and phase transformation stress by FES reveals the origin of residual 
stress distribution feature in quenchedzeylinders, is, the relative higher phase transformation compressive 
stress and lower tensile thermal stress at the core of cylinder makes the residual stress to be compressive, while at 
the surface of cylinde compressi % stress is predominantly from thermal stress, because it is much larger than 
the tensile stress NÉ phase transformation stress. The tangential residual stress distributions in cylinders 
with several diameters 人 人 om 3 mm to 100 mm were predicted by FES, and the results indicate that when diameter 
is less dis mm, the tessile stress at the surface increases with increasing diameter until to 5 mm, then decreases 
with increasing diameter to 20 mm, finally the tensile stress bec compressive stress. Besides, with the 
increase of diameter, maximum tensile stress shifts from the surface tó the location of 0.6 radius. The effects of 
different quenching media on quenching stress woe investi by FES. The results demonstrated that 
although there is compressive stress at surface of cylinder quenched in water or salt solution, the maximum stress 
locates the subsurface, meaning that cr g a s at the subsurface, which is consistent with cracking in 
practical components. This work is M qe adds the analysis of cracking from quenching stress in components with 
different sizes and under different quere media. 
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Zt VE CI reg da [8] JR] FP RC AERRUS DG REI Z GIA TERR. D uM dak 泛 应 用 于 机 械 制 造 行 
业 。 滩 火 一 般 需要 快 冷 以 获得 马 en | 与 心 部 的 温差 较 大 ， 而 且 奥 氏 体 向 马 氏 体 的 转变 体 脱 
厂 余 应 力 娩 开裂 倾向 。 科 研 人 员 估计 有 20% 的 热 
。 因 此 ， 对 于 溢 火 冷 却 过 程 中 的 应 力 及 变形 的 控 
0 康 大 万 等 中 利用 X 射线 应 办 分 析 法 恨 法 研究 了 调 质 大 轴 残 余热 应 力 分 布 规律 和 柱 形 件 
残余 组 织 应 力 ( 相 变 应 力 ) 的 分 布 规律 。 张 海外 使 用 同样 方法 测量 了 直径 60 mm 和 100 mm 的 18Cr2Ni4W 柱 
形 件 中 残余 应 力 沿 截面 的 分 布 。Wang 等 中 通过 x 射线 方法 测量 了 1080 圆 棒 试 样 在 水 和 两 种 有 机 物 泽 水 介 
质 分 别 滩 火 后 的 残余 应 力 截面 分 布 。 尽 管 可 以 通过 实验 方法 测量 简单 状 试 样 的 应 力 分 布 ， Mrs c 
几 个 限制 : 1)9: 358 WU EEREERERE GE, ARRENE, MEWE, LÍERBUKITAYM ACUAA ELTE 
力 分 布 的 测量 较为 困难 ;3) 实 验 测 得 的 简单 小 试 样 的 结果 无 法 直接 应 用 至 实际 大 型 复杂 工件 ，4) 根 据 目 
前 实验 手段 ， 尚 无 法 测量 出 济 火 过 程 的 瞬时 应 力 (过 程 应 力 )， 而 很 多 情 TURROIA ARETE 
却 过 程 中 发 生 的 号 ，5) 通 常情 况 下 ， 无 法 通过 实验 分 别 确定 淳 火 过 程 中 间 时 存在 的 热 应 力 和 组 织 应 力 ， 但 
其 分 离 在 有 限 元 模拟 中 可 以 很 容易 实现 。 比 如 ， 可 以 通过 分别 计 和 ARQUATA, RERE- HFK 
EUIS AGBRUSL I IRL ESD, AMY AAE HERE ECT ERSA, 调整 工件 尺寸 结构 或 材料 成 分 
等 ) 降 低 该 位 置 最 大 应 力 , 避免 工件 开裂 。 基 于 以 上 原因 , 计算 机 模拟 广泛 地 被 应 用 于 湾 火 应 力 的 研究 1。 
通过 计算 机 模拟 ， 可 以 预测 出 任何 形状 及 任意 大 小 工件 的 温度 、 组 织 及 应 力 的 变化 过 程 。 Simsit s 等 9 使 用 
了 有 限 元 方法 模拟 了 三 维 C60 偏心 圆 环 的 温度 场 、 组 织 场 和 应 力 场 的 变化 ， 并 与 实验 结果 进行 了 比较 ， 结 
果 表 明 ， 建 立 的 有 限 元 模型 能 够 预测 出 三 维 非 对 称 工件 残余 应 力 分 布 趋势 。Arimoto a 
方法 分 析 了 偏心 圆 环 件 开 裂 的 原因 ,计算 的 最 大 主 应 力 位 置 与 开裂 位 置 相 同 。Jung 等 中 提出 了 一 种 新 的 确 
定 临 界 温度 变量 ( 贝 氏 体 开始 转变 温度 (Bs) 和 马 氏 体 开始 转变 温度 (Ms) 的 方法 )， 该 改进 的 方法 提高 了 相 变 及 
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应 力 计算 精度 。 Lee 等 1 提出 了 一 种 计 入 
透 性 的 中 碳 调 质 钢 湾 火 应 力 及 组 织 分布 ， 表 面 应力 的 计算 结果 与 XRD Ma 
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Ms 和 马 氏 体 相 变动 力学 的 新 方程 。Ariza ^: 


等 呈 模 拟 了 三 种 不 同 深 
上 量 结果 吻合 较 好 。 


TEKMI 


计算 机 模拟 中 ， 仍 有 很 多 问题 尚未 得 到 较 好 解决 ， 其 中 相 变 塑 性 是 阻碍 热处理 模拟 技术 发 展 的 较为 棘手 的 


难点 之 一 。 相 变 塑 性 是 在 应 力 及 相 变 同时 


于 母 相 屈服 强度 也 会 产生 相 变 塑性 。 在 

本 课题 组 在 Taleb08 的 实验 测定 的 术 
归 一 化 函数 (exponent-modified (Ex-Modified) oe 
函数 中 更 好 地 拟 合 了 Taleb 的 实验 曲线 , 以 致 更 精 古 uda 了 两 个 不 同 直径 42CrMo 圆柱 件 
经 过 40 多 年 的 发 展 ， 尽 管 热 处 理 
理 模拟 精度 的 因素 如 相 变 塑 性 己 汪 、 


> + >y [15-17 
成 共识 ! i, 


待 进一步 提高 。 
在 生产 中 


淳 火 裂 纹 的 


an 


ioi FUSLZDSET ROS TETEUE AUS 73 BUE RS 
组 织 场 和 应 力 场 相互 耦合 的 数学 模型 ， 然 后 通过 测量 


ee 最 后 通 


Ea 度 场 、 


余 应 力 ， 对 模型 的 可 靠 1 性 (包括 相 变 塑性 
因素 进行 分 析 。 


1 实验 方法 


实验 所 用 材料 为 40CrNiMo， 其 
圆 棒 试 样 直径 为 60 mm/s 
首先 将 试 样 放 入 温度 为 850 
出 ， 预 冷 30 s 后 迅速 垂直 放 入 水 温 
中 一 个 用 来 进行 组 织 及 硬度 分 析 ， 


k, 


发 点 25 


使 用 线 切 割 将 其 中 一 
行 抛光 ， 抛 光 后 将 样 
CeH3N307+2.5 ml HCI+50 ml C 


NAE rp 


rH m 
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技术 在 


善 ， 并 且 已 7 
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另 一 个 用 于 应 力 测量 。 
面 切 开 ， 在 1/2 半径 位 置 
JÆ 2%(2 ml HNO4*98 ml C 
中 腐蚀 10 


PHERI P, 


EH 


NHE 


保温 100 min 使 样品 温 
为 室温 (32 c) kir, AA 150 s ERETTA HAH 


a UB 


吹 干 试 样 。 将 处 理 好 的 样品 放置 于 金 相 显微镜 (OM， 府 司 ， 


软件 中 的 快速 选择 工具 
所 占 像素 量 除 以 整 引 


勾 选 出 金 相 图 片 
KE ER SR ECRIRE T 


v 的 少量 


A. AKKA 

Imager Alm) 中 进行 观察 。 
EI Pin 量 贝 氏 体 组织 并 计算 其 所 占 像素 量 
n FPE, RRAS REHAR. Vickers 硬度 EE 


GED e HEX EISE 


青水 和 酒精 


存在 的 情况 下 产生 的 可 被 观察 到 的 塑性 应 变 ， 即 使 


其 等 效应 力 低 


述 工件 济 火 过 程 中 的 应 力 及 应 变 状态 时 ， 相 变 塑 怕 

塑性 的 基础 上 ， 提 出 了 一 种 指数 
ed function)， 比 Abrassarts 和 Desalo 提出 的 相 变 塑 性 
KRKE, 
FREE FEAA, 
点 力 对 相 变 的 影 2g 目前 尚未 得 到 较 好 的 解决 ， 热 处 理 模 拟 技术 仍 


多 式 修正 的 相 变 塑 性 


应 该 被 考虑 已 


由 于 影响 热 处 


， 纵 A NUN 常 Tr RESHKE REREN IMA YEARUE SS, HARKEN 
e a eM E 


IEZA IEEE E 


究 较 少 PC 。 


见 被 其 他 学 者 报导 。 本 工作 中 综合 
Moab ERU 
应 力 分 布 的 特点 及 影响 


为 Fe-0.38C-0.27Si-0.62Mn-0.81Cr-1.35Ni-0.18Mo( 质 量 分 数 ，%)。 


匀 一 致 ， 然 后 将 其 取 


HVS-30P 硬度 仪 下 进行 ， 加 载 载荷 30 kg， 保 持 时 间 15 s， 每 个 位 置 重 复 测 量 3 次 。 


土 37” 和 土 43° 


式 中 el 为 与 试 样 表面 法 线 之 间 的 夹 角 为 wy 的 (211) 衍 射 面 的 法 向 应 变 ; 1/2 S 25.9210 MPa". 


每 个 半 和 en N 


X 射线 仅 能 


量 试 样 表面 残余 应 力 ， 若 测试 样 沿 


使 用 这 RD combo X 射线 应 力 分 析 仪 (Proto 公司 , 加 拿 大 ) 对 试 样 进 4 
为 试 样 表面 法 线 与 衍 
HE Ka 辐射 0=0.2291 nm), 》 
。 残 余 应 力 的 计算 公式 为 : 


VR 


射 晶 面 法 线 之 间 的 夹 角 )， 测 量 执行 EN15305-2008 标准 27。 使 
ur MM 


{211} 
0g, 
Osin?j 
QU 


h 


REANIM, 


了 应 力 测量 , 测量 方法 为 sip2y 法 (w 
EL oM. 


F 经 粗 磨 精 磨 后 进 


E 2 个 试 样 ， 其 


HRW RONA Vilela 试剂 (0.5 g 
清洗 ， 最 后 用 冷风 
首先 使 用 Photoshop 


; 贝 氏 体 


«12 | +24° 、 士 30* 


(1) 


腐蚀 后 


须 对 试 样 进行 剥 层 。 为 了 提高 


剥 层 效 
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率 ， 先 使 用 车 床 车 去 0.5 mm 深度 ， 然 后 再 用 化 学 腐蚀 方法 去 除 车 床 加 工 引 起 的 附加 应 力 层 ， 为 完全 去 除 
附加 应 力 层 ， 每 次 腐蚀 的 深度 大 于 0.5 mm。 所 使 用 的 化 学 腐蚀 剂 配 比 为 : 150 ml HNO3+50 ml HzO2+20 g 
C2H2O04， 加 水 至 500 mle5。 剥 层 后 ， 试 样 残余 应 力 会 重新 分 布 ， 为 了 考虑 去 除 材 料 部 分 对 所 测 位 置 残余 
应 力 的 影响 ， 使 用 下 列 公 式 对 其 进行 修正 3; 


R d 
«(r)-e. [^ (0). 52 


» R dé 
ce,(r)-o,.(r)-2]| Orm é) -一 
goo os es 


式 中 ，a(D 和 oi(n) 分 别 为 半径 为 r 处 修正 后 的 切 向 应 力 和 轴 向 应 力 ，oin() 和 oon(n) 分 别 为 修正 前 的 切 向 应 
JARIR: RIPEKA m) EAEE RAAE ERE 
2 woe 

str, AN ERAN, SED ARRIERE RI tto, R 
管 应 力 对 相 变 动力 学 有 影响 64, [8 erbe RT AE DU S HEX RI REG 在 本 工作 中 暂 不 考虑 该 交互 作用 。 
文献 9 指出 在 淳 火 过 程 中 因 工 件 变形 而 引起 的 温度 变化 不 超过 2 它 ， 因 此 , 为 了 简化 ,该 影响 在 本 工作 中 
亦 不 考虑 。 第 三 个 简化 是 假设 材料 成 分 均匀 。 根 据 水 比 简 化 ， 泊 | 医 刀 、 组 织 场 和 应 力 场 之 间 的 全 厅 合 关系 
可 简化 为 顺序 耦合 关系 (如 图 1 所 示 )。 


(em 
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Fig.1 Simplified oE of quenching process 
5 E 
24 温度 场 计算 P 


因 不 考虑 塑性 变形 对 温度 的 影响 ， 整 个 温度 场 计算 的 支配 方程 可 写成 9 5 
pC, (o, T)T 2 V(k(o,, T). VT)-Q UP K (3) 


式 中 ，p、Cb 和 K 分 别 为 混合 组 织 的 密度 、 比 热 容 和 热 导 率 ; 了 为 温度 ; 元 为 温度 对 时 间 的 偏 导 ; YY 为 梯度 
算 子 ;0@ 为 单位 时 间 内 相 变 潜 热 释放 量 ， 其 表达 式 为 : € 


à-Ynfgy, ES (4 


式 中 , AH, IEK; px 为 k 组 织 (k=1, 2, 3 分 别 代 表 奥 氏 体 、 贝 氏 体 和 马 氏 体 ) 的 体积 分 数 ; 73 k 2H ZR CUL 
氏 体 或 马 氏 体 ) 的 相 变 速率 。 严 格 来 说 , p 是 组 织 和 温度 的 函数 ， 因 其 微小 变化 对 温度 场 计算 的 结果 影响 较 
小 ， 因 此 在 温度 场 计算 过 程 中 jp 取 一 恒定 值 : 7850 kg/m’. C, k 的 值 取 自 文献 5 提供 的 数据 。 奥 氏 体 向 
贝 氏 体 或 马 氏 体 转变 的 相 变 答 分 别 为 -5.12x108J/m2 或 -3.14x108J/ms03。 

工件 的 初始 温度 设置 为 850C 。 工 件 表面 换 热 的 边界 条 件 为 : 


(2) 
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-k(puT)VT =h(T HT, -1,) 6) 


式 中 ， ee 其 值 是 工件 表面 温度 的 函数 ， 五 和 dod eui ce RE it E o 

于 影响 水 的 换 热 系数 的 因素 较 多 (如 水 温 、 水 槽 大小、 搅拌 速率 和 材料 表面 状态 等 )， 因 此 不 同文 献 
Eu 00 为 了 降低 因 温 度 计算 误差 对 应 力 计 算 结 果 的 影响 ， “本 工作 对 直径 60 
mm 的 42CrMo 圆 棒 离 表面 5 mm 位 置 、1/2 半径 位 置 和 中 心 位置 的 水 冷 温度 进行 了 测量 。 根 据 所 测 得 的 冷 
却 曲 线 拟 合 ， 对 水 的 换 热 系数 进行 了 估算 。 WC Him 具体 步骤 如 下 : 1 将 换 热 系数 离散 为 
干 温度 点 ， 两 个 相 邻 温度 点 之 间 的 换 热 系数 采用 线性 内 插 法 计算 ; 2) 对 这 些 离 散 温度 点 设 初 始 值 ， 其 初始 
值 可 从 相关 文献 中 查 到 ”2,; 3) 采 用 粒 乞 群 优化 方 夸 lartiaqle swarm optimizing method) 对 这 些 温度 点 的 换 热 系 
数值 进行 优化 ， 在 该 方法 中 目标 函数 为 计算 得 到 的 冷却 曲线 与 实测 冷却 曲线 围 成 的 面积 ，3) 根 据 目 标 函 数 
修改 换 热 系数 ,直至 escis 的 换 热 系数 及 计算 冷却 曲线 与 实测 冷却 曲线 的 比较 如 图 2 所 示 。 


En 


900 14000 
( sop” —— Core(Sim.) G 
> S + Core(Meas.) 人 12000 | 
| ---12R(SimJ) [S 
12R(Meas ) & 10000 - —m— Modified 


—-—— 5mm from surface(Sim.) —e-— Ref. 29 
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E 400 22 E 6000 - 
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E z 
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200. E 
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图 2 "TTL TT T 结果 与 m ARE 
数 的 比较 


Fig.2 Comparison of measured and calculated cooling curves (a), and comparison of 2. Nd heat transfer 
coefficient and that from Ref. "2 


22 ”组 织 场 计算 X 


40CrNiMo JB PURGE, XE 60 mm "E 并 无 铁 素 体 或 珠光 体 转变 ， 因 此 本 工作 
仅 考虑 由 氏 体 和 马 民 体 生变 - SEP EAE, It KEISERE, 通过 K-M 公式 9 对 其 进行 计算 ， 


ov -1-exp| -a(M, -T;) | XS. 6 
RH, gu 为 马 氏 体 体积 分 数 ，M4 为 马 氏 体 开始 转变 温度 ， 0 qe 
时 最 低 冷 却 温度 ，4 为 一 系数 ， 其 数值 取决 于 钢 中 化 学 成 分 3 


à: —0.0224-0.0107x.. — 0.0007 Xm 一 0.00005xNi — 0.00012x.., — 0. qp. (7) 
贝 氏 体 相 变 具有 扩散 型 相 变 特征 ， 在 恒温 情况 下 ， 其 相 变 量 通过 CONES 


Pr -1-00 $ l 3$ (8) 


AP, pe 为 贝 开 体 体积 分 数 ; zf 为 转变 63.2% 贝 民 体 所 需 时 间 ; n 为 指数 因子 。 为 了 将 公式 (8) 应 用 在 连续 
冷却 情况 下 ， 使 用 Scheil 姜 加 法 则 中 的 方法 进行 计算 ， 其 具体 计算 过 程 可 参考 文献 ，。 


2.3 ”应 力 场 /应 变 场 计算 
因 本 工作 中 所 研究 的 工件 尺寸 较 小 ， 滩 火 持 续 时 间 短 ， 在 此 不 考虑 蠕 变 影响 ， 故 而 工件 中 任意 位 置 的 
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总 应 变 增 量 Ae; 可 分 解 为 : 


Ag, = Ae? t Ag? 4 Ag! HAE" AZ? (9) 
ij ij ij ij y 


RP, Aeh, Ae, Ae”, Ae" 和 Ae" 分 别 为 弹性 应 变 增 量 、 塑 性 应 变 增 量 、 热 应 变 增 量 、 相 变 应 变 增 
量 和 相 变 塑性 应 变 增 量 。 A 
热 应 变 和 组 织 应 变 属 于 各 向 同性 应 变 ， 其 计算 这 程 如 下 


29 m *AT)- 4 (T)] 
k=1 (10) 


€ Nr he [st (T - AT)- & (T - AT)] 


AF, KEK 2, 5 ca 贝 氏 体 和 马 氏 体 ，gx 为 k 组 织 的 相对 应 变 。 假 设 奥 氏 体 在 0 "C 时 的 热 应 
变 为 0， 根据 文献 B9 所 提供 的 热膨胀 系数 及 我 们 所 测量 的 膨胀 CI T 计算 出 不 同 温度 下 奥 氏 体 、 贝 氏 体 
和 马 氏 体 各 自 的 膨胀 曲线 ， 如 图 3 所 示 。 

cq 


oook « 
pe ——— Austenite E 


- - - Bainite 


0.005 - 一 一 一 B^ K 
0.000 1 1 L L L L L K 
O 100 200 300 400 500 600 700 800 900 


"C 


图 3 奥 氏 体 、 贝 氏 体 和 马 氏 体 的 热膨胀 曲线 

Fig.3 Thermal expansion curv ustenite, bainite and martensite 
当 相 变 在 应 力 (即使 低 于 母 相 屈服 强 co 下 发 生 时 ,材料 会 产生 异常 塑性 , 该 塑性 被 称 ME ppl, 
相 变 塑性 增 量 的 表达 式 采用 本 课题 组 提出 的 相 变 塑性 函数 中 


t 3 n n= 
^c? => AmnkK (c exp(mg e '"AQS; zz Loo (11) 


AF, A. m 和 均 为 拟 合 常数 ， 其 数值 分 别 为 1.105、-0.0994 和 -0.91; K(o) 为 相 变 塑 性 系数 ， 其 值 通 
拟 合 文献 [40] 中 实验 数据 得 到 ;si 为 偏 应 力 张 量 。 €& 
根据 应 变 增 量 ， 可 以 求 出 应 力 增 量 


2G 28. ，, 2A Ee : 2 
nie e] nS o 


p4 =Š Amnk (o)exp(mo")g""Ap (12) 
l=1+2Gy 


Am, jc dud B 为 体积 模 量 ，G 为 剪 切 模 量 ， x 和 1 均 为 变量 。 
混合 组 织 的 屈服 函数 中 关 
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(13) 


AF, oi 是 加 载 曲面 的 中 心 在 应 力 空间 内 的 移动 张 量 。 当 @>0 时 , 材料 产生 届 服 , 使 用 Euler 回 退 算法 [41] 
对 其 进行 塑性 修正 。 混 合 组 织 的 屈服 强度 表达 式 为 : 


e,-8g (9 )o,a 十 | 1- 9 (9, )| O y BM (14) 
n 
式 中 , oy 是 混合 组 织 的 届 服 强度 , ovA 为 奥 氏 体 的 届 服 强度 ; miaw 为 贝 氏 体 和 马 氏 体 的 平均 届 服 强度 ; goa) 
为 奥 氏 体 体积 分 数 的 归 一 化 函数 ， ERNOUR 1142]. 
24 ”模拟 过 程 


本 工作 中 "m N ee 求解 器 进行 求解 , 组 织 场 通 过 修改 子 程序 进行 加 入 。 相 
变 潜 热 在 子 程序 U T 中 加 入 ; 相 变 及 相 变 塑性 在 子 程序 UMAT 中 加 入 。 因 本 工作 所 涉及 的 样品 形状 
均 为 柱 奖 纪 游 火 时 蕊 是 迅速 垂直 入 水 ， 因 此 其 满足 轴 对 称 条 件 ， 故 而 为 了 提高 计算 效率 ， 温 度 场 和 应 力 场 
计算 单元 分 别 选择 DCAX4 和 CAX4 轴 对 称 单元 。 E T ORR EAE, 划分 网 格 时 对 表面 附 
近 的 网 格 进行 了 细 化 。 


3 ”实验 和 模拟 结果 K 
3.1 组 织 和 硬度 分 布 2 


4 为 40CrNiMo 直径 为 60 ad 1/2 半径 和 中 心 位置 的 组 织 OM 像 。 从 图 中 可 
以 看 出 ， 在 心 部 及 12 HEMER RE, EREDA DERN RER: 心 部 贝 氏 体 体积 分 
数 约 8%; 1/2 hr LAE BS TERCER ABUS. 2%。 心 部 位 置 贝 氏 体 含量 高 于 1/2 半径 位 置 的 原因 是 心 部 位 置 
的 平均 冷 速 ( 约 -13 'C/s) 低 于 1/2 半径 位 置 的 平均 冷 速 ( 约 -15CAS)。 因 心 部 贝 eco. 且 该 位 置 维 氏 
硬度 与 表面 位 置 维 氏 硬度 相差 其 小 ， 因 此 可 以 认为 该 试 样 被 完全 淳 透 。 图 5 中 计算 后 的 组 织 分 布 与 实测 结 
果 具 有 相同 趋势 。 图 5 中 硬度 的 分 布 结果 表明 ， iB EHE JR UCET RR 部 K 心 部 位 置 附近 硬度 低 于 
表面 附近 硬度 的 原因 是 心 部 附近 含有 相对 多 的 贝 氏 体 组 织 。 


QNA 


4 60 mm 直径 的 40CrNiMo EEK% 1/2 d c OM 像 


Fig.4 OM images of the quenched 60 mm diameter cylindeVaihe locati: 1/2R (a) and core (b) respectively 
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| B 

m ee ou J 
9—"—9 

C ———— — 


* Bainite(Exp.) 

-- -ABainite(Sim.) 
ensite(Exp.) 
ensite(Sim.) 


—e— HV(Exp.) 


F 水 滩 后 组 织 体积 分 数 ( 计 算 和 测量 ) 和 硬度 ( 测 


量 ) 的 截面 分 布 


图 5 No E 
Fig 5 parison bétween the calculated and measured microstructure fraction distribution, and measured HV 


3.2 ”残余 应 力 分 布 


distribution along radius of octo 


ecu di 
为 了 验证 模型 的 适用 性 ， 比 较 了 直径 60 mm if SOCINIMo BLUES VEEZK IRE SERIAL. W IRRI BS 


XRD 测量 与 计算 结果 ， 如 图 6 Har. WEF, 
同 的 应 力 分 布 特点 : 1) 切 向 最 高 残余 
均 处 于 工件 的 次 表层 位 置 ， 并 非 处 于 玉 件 表面 位 置 ，3) 工 伯 


ge 应力 的 计算 结果 与 模拟 结果 符合 较 好 ， 两 者 均 上 共有 相 
力 高 于 轴 何 最 高 残余 拉 应 力 ; 2) 轴 向 和 切 向 残余 拉 应 力 的 最 大 值 
表面 位 置 为 残余 压 应 力 ， 且 切 向 残余 压 应 力 的 


绝对 值 高 于 轴 向 残余 压 麻 罚 的 绝对 作 \ 为 在 心 部 位 置 ， 工 件 的 轴 向 及 切 向 残余 应 力 均 为 压 应 力 。 应 力 分 布 
与 文献 [426,43,44] 所 提供 的 全 济 透 件 应 力 分 布 特征 相同 。 


4 分 析 讨 论 


-800 - Tangential Stress(Sim.) 
-1000 L 
1200 1 1 1 1 1 1 
0 5 10 15 20 25 30 


Em Tangential Stress(Exp.) 


Radius / mm 


图 6 60 mm 直径 的 40CrNiMo AEJK Ja t [8] 717] 


Fig.6 Axial and tangential stress distribution curves of 60 mm diameter 4 


& Y* 


4.4 直径 60mm 的 40CrNiMo 圆 棒 水 淳 后 应 力 分 布 特征 分 析 


从 图 6 中 残余 应 力 的 分 布 可 以 看 出 ， 切 向 残余 应 力 与 轴 
残余 拉 应 力 高 于 轴 向 最 高 残余 拉 应 力 ， 因 此 全 淳 透 柱状 工 


的 截 


rNiMo bar after water quenching 


S- 


向 残余 应 力 的 分 布 特征 基本 相同 ， 但 切 向 最 高 
人 易于 产生 纵向 裂纹 号 。 故 而 对 于 全 淳 透 性 件 ， 


其 切 向 应 力 的 讨论 较 轴 向 应 力 更 具有 实际 意义 。 为 了 找 出 造成 最 大 残余 切 向 应 力 的 原因 ， 通 过 模拟 绘 出 了 


最 大 残余 切 向 应 力 位 置 ( 距 表 


Ej 10 mm 处 ) 表面 和 心 部 位 置 切 向 应 力 随 冷 却 时 间 的 变化 过 程 , 如 图 7 所 示 。 


在 空冷 阶段 的 30 s A, 表面 为 拉 应 力 ， 心 部 为 压 应 力 ， 这 是 由 于 空冷 阶段 表面 温度 的 降低 引起 的 收缩 使 得 
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工件 表面 为 拉 应 力 ， 由 于 空冷 较为 缓和 ， 工 件 截面 温差 小 ， 此 时 的 应 力 较 小 。 当 工件 入 水 瞬间 ， 因 表面 温 
度 的 迅速 降低 使 得 表面 拉 应 力 迅速 增加 ， 而 心 部 的 压 应 力 增 加 较为 平缓 ， 这 是 因为 心 部 温度 较 高 ， 奥 氏 体 
在 高 温 状 态 下 届 服 强度 较 低 。 当 工件 表面 温度 降低 至 Ms 温度 以 下 时 ， 由 于 表面 马 氏 体 相 变 产生 膨胀 ， 且 
其 膨胀 受到 心 部 的 约束 , 造成 表面 的 拉 应 力 反 向 变 为 压 应 力 , 而 心 部 及 距 表面 10 mm 位 置 因 表 面 的 马 氏 体 
膨胀 而 处 于 拉 应 力 状 态 。 当 距 表 面 10 mm 位 置 处 的 温度 降低 至 Ms 时 ， 其 拉 应 力 降 低 的 同时 引起 心 部 拉 应 
力 的 增加 及 表面 压 应 力 的 降低 。 继 续 冷却 ， 当 心 部 位 置 的 温度 降低 至 Ms 温度 以 下 时 ， 心 部 拉 应 力 降低 逐 
步 转变 为 压 应 力 ， S EEA Oma C TAPER EE SR. Se EE 
变 结束 后 继续 冷却 时 ， 因 心 部 温度 的 继续 降低 ,引起 表面 压 应 力 的 增加 及 距 表 面 10 mm 位 置 处 拉 应 力 的 降 
低 。 从 图 7 中 还 可 以 看 出 ， 下 表面 Bgam 位 置 及 大 切 向 拉 应 力 所 对 应 的 时 间 是 心 部 马 氏 体 完 全 转变 的 


时 间 。 
A 400 } 10 mm from surface 
200 } 
原 人 
z 
j -200- 
à -400 - 
S -600- 
Is 
& -800 上 
[S 
}= -1000 - 
-1200 - 
-1400 - 


1 L L 1 1 L L 1 1 1 L 
0 15 30 45 60 75 90 105 120 135 150 165 180 


tis 
g Ww 
图 7 直径 60mm 的 40CrNiMo AEKA O — ETE I EB e HA 10mm 位 置 处 向 应 力 变 化 过 程 


Fig 7 Evolution of tangential stress with time at the location of core, surface and 10mm below surface of 60 mm 


diameter cylinder of A0CrNiMo during quenching in water 
由 于 汉 火 残余 应 力 是 热 应 力 和 组 织 应 力 两 者 相互 个 加 的 结果 ， 分 别 计算 了 工件 的 热 应 力 、 组 织 应 
力 ， 如 图 8 所 示 。 热 应 力 截面 分 布 的 特点 是 : 表面 压 应 力 ， 心 部 拉 应 力 ; 组 织 应 力 截 分 布 的 特点 是 ; 心 
部 压 应 力 ， 表 面 拉 应 力 。 最 终 残 余 应 力 的 分 布 特点 取 集 过 这 两 者 应 应 的 相对 大 小 。 从 图 8 中 可 以 很 清楚 地 
看 到 工件 表面 及 心 部 的 压 应 力 分 别 是 由 热 应 力 和 组 织 应 力 引起 ,而 距离 表面 10 mm 处 最 大 切 向 拉 应 力主 要 


是 由 组 织 应 力 引起 。 [3 ES 
S- 


900 - 

s 1 
$ 300. K 
= 0 
9 | ~ 
E -300 - K 
[5] 

E: -600 - 
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图 8 直径 60 mm 的 AOCrNiMo 圆 棒 水 济 热 应 力 、 组 织 应 力 及 残余 应 力 计算 结果 的 比较 
Fig.8 Thermal stress, transformation stress and residual stress of 60 mm diameter cylinder of 40CrNiMo after 
water quenching (transformation plasticity was included during all the calculations) 
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K TARARE KARRA, 通过 软件 计算 直径 3 mm. 5 mm. 10 mm, 20 mm, 40 mm, 
70 mm 和 100 mm I 1155 PE ARE ZK UE BRE T8056 28 II 73 V 1 3 7E] 89] 2 2 )WI 7] 438 , 如 图 9 所 示 。 
从 图 9(a) 可 以 看 出 以 下 几 个 规律 : 1 在 直径 小 于 5 mm 时 ， 随 着 直径 的 增加 ， 表 面 残余 拉 应 力 逐 渐 增 加 ， 


当 直 径 大 于 5 mm 时 , 表面 残余 拉 应 力 开始 降低 , 当 直 径 增 加 至 20 mm 时 , 表面 残余 应 力 已 转变 为 压 应 力 ， 


Y 置 3) 截面 最 大 拉 应 力 数 值 随 直 径 的 增加 具 


EE x 


。 图 9(b) 和 图 9(O 表 明 ， 热 应 力 和 组 织 应 力 均 随 着 工件 直径 的 增加 而 增加 ， 但 因 工件 热 应 力 随 直径 的 增 


应 力 数值 随 直径 的 增加 而 逐渐 增加 ;2) 随 着 渤 件 直径 的 增加 ， 最 大 拉 应 力 从 表面 逐渐 过 渡 至 0.6 倍 半 


增加 (<5 mm) 后 减 (5 mm 至 20 mm) 最 后 逐渐 增加 的 趋 


加 增幅 高 于 组 织 应力 ， 因 此 ， 随 着 千 荐 直径 的 增加 ， 表 面 压 应 力 值 逐渐 增加 。 在 0.6 半径 位 置 附近 的 热 应 


力 逐 渐 过 渡 至 0.6 半 


o 


fice cdi p 并 且 该 数值 均 随 着 直径 的 增加 而 增 大 ， 因 此 ， 随 着 工件 直径 的 增加 ， 最 大 拉 应 
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图 9 [8] ÉL Ae VER DL T PEZ AUS PROT E AISISE TT 热 应 力 和 组 给 力 分 布 K 


Fig.9 Distributions of residual tangential stress(a), thermal stress(b), phase tzansformatioivstress(c) of the fully 


harden cylinders with various diameters after quen in water 


4.3 EKNER RSS IRURE LT L7] 43 A6 BERE 


S$ 
l 
选择 4 种 不 同 的 冷却 介质 进行 计算 对 比 研究 ， 人 10%NaCl 水 溶液 、 


EEEE. 10 为 直径 50 mm 圆 棒 试 样 分 别 在 1 
切 向 残余 应 力 、 热 和 组 织 应 力 模拟 结果 。 可 以 看 出 ， 


清水 、80C 淳 火 油 和 空气 。 除 水 的 换 热 系数 取 自 图 2(b) 外 ， 其 他 三 者 的 换 热 系数 取 自 文献 [1]。 假 设 材料 能 


0%NaCl 水 溶液 、 清 水 、80'C 泽 火 油 和 空气 中 汉 火 后 的 
汉 火 介质 对 滩 火 残余 应 力 的 影响 极 大 : AS IRURE AT 


质 济 火 后 的 应 力 分 布 及 最 大 拉 应 力 大 小 及 位 置 完 全 不 同 。 空 冷 后 的 残余 应 力 最 小 , 在 10MPa 以 内 , 这 是 由 


于 空冷 较为 组 和， 热 应 力 和 组 织 应 力 均 较 小 的 缘故 ( 
的 温差 较 小 ， 因 此 因 温 度 梯 度 造 成 的 热 应 力也 较 小 ( 


图 10b)。 对 于 油 冷 ， 由 于 冷却 较为 缓慢 ， 表 面 与 心 部 
图 10b)， 又 由 于 马 氏 体 相 变 膨 胀 引起 的 组 织 应 力 高 于 
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热 应 力 ( 图 10(b))， 因 此 其 合成 应 力 具 有 组 织 应 力 分 布 特点 (图 10a) 。 模 拟 结果 与 文献 [45] 所 提供 的 实测 结果 
相符 。 清 水 和 10%NaCl 水 浴 液 因 冷 速 较 高 ， 表 里 温差 较 大 ， 热 应 力 和 组 织 应 力 均 较 高 ， 淳 火 后 均 具 有 较 
高 的 残余 应 力 ， 其 中 因 10%NaCl 水 溶液 冷 速 高 于 清水 的 冷 速 ， 因 此 其 表面 具有 更 高 的 残余 压 应 力 。 结 合 
图 10a 和 b 可 以 看 出 最 大 残余 拉 应 力 位 于 16 mm 半径 位 置 处 , 该 位 置 的 热 应力 和 组 织 应 力 均 为 拉 应 力 ， 因 
10%NaCl 水 溶液 热 应 力 较 高 ,， 因 此 其 全 加 后 的 总 的 拉 应 力 蜂 值 高 于 水 济 后 的 拉 应 力 峰 值 。 AEKA K 
淳 火 后 其 表面 为 压 应 力 ， 但 其 最 大 拉 应 b 应 力 ， 这 与 文献 [26] 提 供 的 水 淳 裂纹 经 常 发 生 
在 表面 为 压 应 力 的 工件 的 次 表层 的 实验 现象 相符 
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图 10 直径 


Fig. 10 Tangential residual stress 
50 mm diameter fully hardening TAS after quenching in air, 80 °C oil, water and 1096 NaCl solution 


5 结论 S 
(1) ARERR JW 77-5; 40CrNiMo  ZIGAE SUPE IE 4s 8E 73 23 18 B6 UI E vitibus, 证 明了 本 


研究 运用 的 温度 场 、 组 织 场 和 应 力 场 相互 而 合 的 数学 模型 (包括 建立 的 相 变 塑性 函数 ) 的 正确 性 ; 
(2) 对 于 直径 60 mm 的 40CrNiMo RE, Ki sumo E A. 最 大 拉 应 力 位 于 其 次 表层 ， 


且 切 向 拉 应 力 峰值 高 于 轴 向 拉 应 力 峰值 ， 其 最 大 拉 应 力 是 由 工作 心 者 

(3) 根据 工件 直径 的 不 同 ， 淳 透 件 水 淳 后 最 大 拉 感 禾 分 布 在 工 
增加 ， 拉 应 力 峰值 逐渐 过 渡 至 0.6 半径 位 置 ， 

(9 35 50m NEREA EAER 基数 值 小 于 水 和 盐水 济 火 后 的 应 力 ， 通 过 增加 
淳 火 介质 的 冷却 强度 ， 可 使 工件 表面 的 感 为 状态 逐渐 BRUHRBHRERARUS, MEANEN, : 
表面 压 应 力 越 高 ， 但 在 次 表层 容易 引起 开 袭 的 拉 应 力也 越 高 。 

O 对 于 大 直径 工件 、 冷 却 强度 高 的 热 应 力 状态 下 ， 最 高 拉 应 力 位 置 处 于 0.6 半径 位 置 ， 由 于 该 位 置 
热 应 力 和 组 织 应 力 均 为 拉 应 力 ， 因 此 ， 在 该 情况 下 不 能 借助 热 应 力 抵消 组 织 应 力 ， 坷 评 力 和 组 织 人 为 均 是 
促进 工件 开裂 的 原因 。 


马 氏 体 相 变 膨胀 引起 的 ; 
层 或 次 表层 位 置 ， 随 着 工件 直径 的 
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